Nodeless vs nodal order parameters in LiFeAs and LiFeP superconductors 
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High-precision measurements of magnetic penetration depth A in clean single crystals of LiFeAs 
and LiFeP superconductors reveal contrasting low-energy quasiparticle excitations. In LiFeAs the 
low-temperature A(T) shows a flat dependence indicative of a fully gapped state, which is consistent 
with previous studies. In contrast, LiFeP exhibits a T- linear dependence of superfluid density oc A~^, 
indicating a nodal superconducting order parameter. A systematic comparison of quasiparticle 
excitations in the 1111, 122, and 111 families of iron-pnictide superconductors implies that the 
nodal state is induced when the pnictogen height from the iron plane decreases below a threshold 
value of ^ 1.33 A. 



There is growing evidence that the superconducting 
gap structure is not universal in the iron-based supercon- 
ductors [l| . In certain materials such as optimally doped 
(Ba,K)Fe2As2 and Ba(Fe,Co)2As2, strong evidence for 
the fully gapped superconducting state has been observed 
from several low-energy quasiparticle excitation probes 
including magnetic penetration depth [21, Isl] and thermal 
conductivity measurements [4, 5]. In contrast, signif- 
icant excitations at low temperatures due to nodes in 
the energy gap have been detected in several Fe-pnictide 

Derconductors. These include LaFePO {Tc = 6K) 
BaFe2(As^)2 {T^ < 31 K) ^^, and KFe2As2 
(Tc = 4K) ^12|-[l^. It is quite extraordinary that such 
distinct pairing states appear in closely related members 
of the same class of superconductors. To understand the 
mechanism of superconductivity in iron-based supercon- 
ductors, it is essential to identify what determines nodal 
and nodeless states |15l-[2ll|. 

Theories based on antiferromagnetic spin fluctuations 
suggest that the pnictogen height hpn above the iron 
plane (see Fig.lTJa)) is an important factor in determin- 
ing the structure of the superconducting order parame- 
ter [l6|, [18| - [2Q| . Generally, hpn is much shorter for the P 
based iron-pnictides in comparison to their As counter- 
parts, so a good test of the theory would be to system- 
atically compare As and P based superconductors. Al- 
though this can be achieved in part in the BaFe2(As,P)2 
series, the fully As containing end member BaFe2As2 is a 
nonsuperconducting antiferromagnet. The same is true 
for LaFeAsO which is the As analogue of the nodal su- 
perconductor LaFePO. Charge doping of the arsenides 
induces superconductivity, but also introduces disorder 
which complicates the identification of the pairing state. 

The 111 materials, LiFeAs [H, |23J and LiFeP [H, [25J 

provide a unique route to study this problem as both ma- 
terials are superconducting {Tc ~ 17K and 4.5 K, respec- 
tively), nonmagnetic, and importantly very clean, with 
long electronic mean- free paths. In LiFeAs, antiferro- 
magnetic fluctuations have been observed [26, 27] and 
fully gapped superconductivity has been demonstrated 
by several experiments [28l-l33j, but no information has 



been reported for the pairing state in LiFeP. Band- 
structure calculations show that the two materials exhibit 
similar Fermi surface shapes [3J, [35|; quasi-cylindrical 
hole sheets near the zone center and two warped electron 
sheets near the zone corner. 

Here we report on precision measurements of the mag- 
netic penetration depth A(T) in single crystals, which 
demonstrate a nodal gap state in LiFeP in sharp con- 
trast to the nodeless state in LiFeAs. Our analysis based 
on accumulated A(T) data in the 1111, 122, and 111 se- 
ries of superconductors, indicates that the nodal state is 
induced when hpn is below a threshold value. By com- 
paring calculated electronic band structures of LiFeAs 
and LiFeP, we discuss the origin of this behavior. 

Single crystals of LiFePn {Pn = As or P) were grown 
by a flux method [36]. The crystal size of LiFeP is up 
to 135 X 135 X r^ 10 /im, which is smaller than that of 
LiFeAs. To avoid degradation of the sample due to re- 
action with air, the crystals were handled in an argon 
glove box and encapsulated in degassed Apiezon N grease 
before measurements. Large residual resistivity ratios 
(- 50 for LiFeAs and - 150 for LiFeP) [36], observa- 
tions of de Haas- van Alphen (dHvA) oscillations in mag- 
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FIG. 1. (Color online), (a) Schematic crystal structure of 
LiFePn {Pn — As or P). The arrow defines the pnictogen 
height hpn from the iron plane, (b) The ac susceptibility 
of LiFePn crystals measured from the frequency shift of the 
TDO. 
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FIG. 2. (Color online), (a) Low-temperature change in the 
magnetic penetration depth in single crystals of LiFeAs and 
LiFeP. (b) Temperature dependence of normalized superfluid 
density Ps{T). We used A(0) = 210 and 150 nm for LiFeAs 
and LiFeP, respectively, (c) Expanded view of ps (T) of LiFeP 
at low temperatures. The solid line is a fit to the T-hnear 
dependence. 



netic torque [37], and sharp superconducting transitions 
(Fig.[Hb)) show that the crystals are of very high quality. 
The temperature dependence of change in the magnetic 
penetration depth was measured by the tunnel diode os- 
cillator (TDO) technique d, Q down to T/T^ ^ 0.03. 
A weak ac field is applied along the c axis so that the 
super current flows in the ah plane. 

Figure[2l^a) depicts the low-temperature variation of 
the in-plane penetration depth AA(T) = A(T) — A(0). 
The data for LiFeAs is completely flat within the ex- 
perimental error of ~ 0.3 nm below T/Tc ^ 0.1. This 
demonstrates negligible quasiparticle excitations at low 
temperatures, indicating a fully gapped state. This result 
is fully consistent with previous results in LiFeAs [29|, |30| . 
In sharp contrast to this, the data for LiFeP exhibits 
much steeper temperature dependence of AA(T) at low 
temperatures. When we use a power law fit AA(T) ex T" 
to this data below T /Tc ~ 0.25, we obtain a small value 
of Of ^ 1.3. In iron-based superconductors, a power law 
dependence with a ^ 2 can be expected even in the dirty 
full gap case when the sign changing s± state is consid- 
ered 1^] , and indeed a tendency of the exponent decrease 
from Of > 3 to ~ 2 with increased impurity scattering has 
been observed experimentally [2|,l39|. However, the small 
power a < 1.5 cannot be explained by such a dirty node- 
less state, and it is rather a strong indication that the 



superconducting gap has line nodes. Indeed, our data 
can also be fitted to oc T'^ / {T -\- T^) ^ which is applica- 
ble to the nodal case with small impurity scattering [40| . 
The obtained low value of T* ^ 0.3 K indicates a clean 
nodal behavior and is consistent with the other measures 
of sample quality described above. 

We also analyze the normalized superfluid density 
Ps{T) = A2(0)/A2(T) (Fig.IStb)). To do this we need 
the value of A(0), which we cannot directly determine 
from the TDO measurements. The small angle neutron 
scattering measurements of LiFeAs reveal A(0) ~ 210 nm 
[28|. To estimate the A(0) value for LiFeP, we consider 
the difference of the effective mass in these two supercon- 
ductors whose carrier number (Fermi surface volume) is 
quite similar. The effective masses determined by the 
dHvA oscillations [37] as well as the electronic specific 
heat coefficients 7 have a factor of ~ 2 difference (7 ^ 16 
and ^ 30mJ/K2mol for LiFeP [24] and LiFeAs |41|], re- 
spectively), from which we estimate A(0) ~ 150 nm. The 
extracted temperature dependence of Ps{T) shows con- 
trasting behaviors for As and P cases at low temperatures 
again: flat dependence for As and steeper dependence for 
P. The expanded view at low temperatures (Fig.[2fc)) 
demonstrates a wide temperature range of T-linear de- 
pendence, which clearly indicates the energy-linear den- 
sity of state of quasiparticles and hence the existence of 
line nodes in the energy gap. 

The strength of the electron-electron correlations can 
be measured by the mass enhancement which is closely 
related to the 7 value. The larger 7 for LiFeAs than for 
LiFeP suggests weaker correlations in the P case, which 
is reinforced by the smaller A value of the Fermi- li quid 
coefficient in the AT^ dependence of resistivity [36|, |42| 
and smaller quasiparticle mass enhancements measured 
by quantum oscillations [37]. Strong correlations usu- 
ally promote sign change in the superconducting order 
parameter [14], which leads to the gap nodes in single- 
band superconductors. In the present multiband case 
with separated Fermi surface sheets, however, the seem- 
ingly opposite trend that LiFeP has nodes but is weakly 
correlated suggests that other factors are also important 
for node formation. 

To obtain further insights, we gather the available 
data for the low-energy quasiparticle excitations in sev- 
eral iron-pnictide superconductors including 1111 p, 0, 
|43L |44, 122 |2|, 13|, l9|, [ij, |45|, and Ill-based materials 
[29|, [30|. The amount of thermally excited quasiparti- 
cles is directly related to the change in the penetration 
depth. Thus we quantify AA(0.2Tc) as a measure of exci- 
tations so that we avoid ambiguity resulting from uncer- 
tainties in A(0). Among the available data in the litera- 
ture, we select only the data which shows either a < 1.5 
or > 2.5 in the power-law approximation, because the 
power-law dependence with a ^ 2 cannot distinguish the 
dirty nodeless and nodal states as discussed previously. 
A plot of AA(0.2Tc) as a function of pnictogen height 
hpn in Fig. [3] suggests that there is a threshold value of 
~ 1.33 A, below which all the superconductors exhibit 
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FIG. 3. (Color online). Pnictogen-height dependence of 
AA(T) at T = 0.2Tc as a measure of low-temperature quasi- 
particle excitations in 1111 (diamonds), 122 (circles), and 111 
(squares) series of iron-pnictide superconductors. Here the 
pnictogen height hpn (see Fig.[TJa)) is determined from the 
structural analysis at room temperature QQ, \^]5^ . 



significant quasiparticle excitations (with a < 1.5) char- 
acteristic of a nodal state. Above the threshold, most 
of the materials are nodeless with the exception of the 
highly hole-doped compound, KFe2As2. This particular 
material is unusual in that it lacks electron sheets and 
thus there is no interband nesting. In additioiLj;he quasi- 
particle effective mass is strongly enhanced [ij, Egj , and 
Tc is very low (~ 4K) so superconductivity may have a 
different origin to that in the other materials. Therefore, 
our analysis strongly suggests that the pnictogen height 
is an important parameter that determines the gap struc- 
ture in the iron-pnictide superconductors having signifi- 
cant interband scattering. One may also ask about the 
Fe-Pn-Fe bond angle, but the nodal LiFeP has a closer 
angle (108.6°) to the perfect tetrahedron value of 109.47° 
than the nodeless LiFeAs (102.8°), from which we do not 
find any simple correlation between the bond angle and 
gap structure in iron-pnictide superconductors. 

The importance of the pnictogen height hpn on the su- 
perconducting order parameter in iron-pnictides has been 
suggested in theoretical considerations based on the an- 
tiferromagnetic spin fluctuation mechanism [16, 18-20]. 
When hpn is low, one of the hole bands with dxy or- 
bital character, which is located near the (tt, tt) position 
in the unfolded Brillouin zone (BZ), tends to sink below 
the Fermi level. The disappearance of this Fermi sur- 
face makes interband electron-hole scattering weaker and 
hence the importance of the scattering between electron 
sheets relatively greater, promoting a sign change of the 
superconducting gap (and hence nodes) on the electron 
sheets. 

To check this, we have performed band-structure cal- 
culations based on density functional theory (DFT) in- 




FIG. 4. (Color online). Calculated Fermi surfaces of LiFeAs 
(a) and LiFeP (b). The spin-orbit interaction is included in 
the DFT calculations. Color indicates the relative weight of 
the dxy orbital contribution, which has been obtained from 
the Wannier fit by using wannier90 '54\ via wien2wannier 
interface [551]. 



eluding spin-orbit coupling, by using the wien2k package 
[47] and the experimental lattice constants and internal 
positions [23, 36]. The obtained Fermi surfaces (Fig.|4]) 
are similar to previous calculations [35| . Although angle- 
resolved photoemission results have suggested quite dif- 
ferent Fermi surfaces with no interband nesting in LiFeAs 
|33| , more recent bulk measurements of dHvA oscillations 
reveal quasi-nested hole and electron sheets ^] in a good 
agreement with the calculations. Importantly, the calcu- 
lations show that the dxy hole sheet, which in the folded 
BZ is the outermost hole sheet at F, is present in both 
compounds. This indicates that the absence of the dxy 
hole sheet is not a requisite for the nodal state. 

A more detailed comparison between LiFeAs and 
LiFeP reveals that the size of the outer hole sheet shrinks 
and its relative weight of the dxy orbital contribution is 
significantly suppressed for the P case. Moreover, the 
middle hole sheet in LiFeP has rather mixed dxy and 
dzx/yz contributions. Such differences in the orbital char- 
acter in hole sheets may affect relative importance of 
the interband hole-electron scattering compared with the 
scattering between electron bands. dHvA measurements 
in LiFeP [37] imply that the middle hole sheets has less 
mass enhancement than the electron sheets, suggesting 
that scattering between electron sheets is stronger than 
the interband electron-hole scattering. This could lead 
to the formation of line nodes in the electron sheets. We 
note that the extended-5 state with line nodes in the elec- 
tron sheets has been discussed as the most likely nodal 
gap structure of BaFe2(As,P)2 [11]. The strong T-linear 
dependence of superfluid density in LiFeP is consistent 
with the nodes being on electron bands containing high 
Fermi velocity parts, which almost coincide with the dxy- 
dominated regions (yellow parts of the electron sheets in 
Fig.JU^b)) [36]. To determine the exact node locations 
in LiFeP, however, other measurements are necessary in- 
cluding angle-resolved probes of low-energy quasiparticle 
excitations. 

Our results that the nodal state is favored for low hpn 



support the trend that the spin-fluctuation theory pre- 
dicts, but there remains chahenging issues including the 
fact that the emergence of nodes is not directly caused 
by the disappearance of the dxy hole sheet. It has also 
been theoretically suggested that a competition between 
the orbital fluctuations and spin fluctuations generates 
nodes in the electron sheets [21]. The difference in the 
orbital character in hole sheets would also change the or- 
bital fluctuations, which may affect the competition and 
hence the nodal gap structure. Further quantitative cal- 
culations of the pnictogen-height effect based on these 
theories will help clarify the mechanism of iron-based su- 
perconductivity. 

In summary, we have measured the penetration depth 
in clean crystals of LiFeAs and LiFeP. We found a T- 
linear superfluid density for LiFeP indicating a nodal 



order parameter, in strong contrast to the fully-gapped 
state found for LiFeAs. A comparison of low-energy exci- 
tations across the different iron-pnictide superconductors 
suggests that the nodal state is induced when the pnic- 
togen height is shortened below a threshold value. 

After completion of this study, we become aware 
of recent thermal conductivity results [56] which sug- 
gest a nodal state in Ba(Feo. 64^110. 36)2As2 with hpn ~ 
1.317 A[57], which supports our conclusion. 
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